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ABSTRACT 

The  forces  on  a  rotor  due  to  asymmetric  pressure  distributions  resulting  from  a 
single  gland  non-circular  labyrinth  seal  in  a  circular  outer  casing  are  analyzed  for  the 
purpose  of  understanding  the  possible  causes  of  synchronous  vibration  due  to  seal 
intolerance.  The  sealing  knife  imperfections  causing  the  non-circularity  may  be  due  to 
manufacturing  defects  or  in  service  ware.  A  lumped  parameter  model  is  developed  for  flow 
in  the  azimuthal  direction  inside  the  seal  gland.  The  resulting  continuity  and  momentum 
equations  are  solved  using  a  regular  linear  perturbation  technique.  Results  from  this  model 
indicate  under  what  conditions  seal  imperfections  can  generate  forces  of  the  same  order  of 
magnitude  as  those  arising  from  residual  rotor  mass  unbalance. 
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I.        INTRODUCTION 

The  United  States  Navy  has  extensive  experience  with  gas  turbine  engines  in 
aviation  applications,  but  only  relatively  recently  has  this  engine  type  been  expanded  to  the 
marine  environment.  Part  of  the  appeal  of  this  propulsion  system  is  its  quick  reaction  time, 
reliability,  maintainability,  and  power  to  weight  ratio.  To  be  successful  on  board  a  Naval 
combatant  the  power  output  must  be  high,  resulting  in  an  increase  in  machine  size  and 
weight  compared  to  typical  aviation  applications.  In  marine  applications  gas  turbine  engines 
are  usually  coupled  in  some  way  to  a  large  set  of  reduction  gears,  that  not  only  reduce  the 
high  rotational  speed  of  the  turbine  down  to  an  applicable  level  for  shaft/propeller  type 
propulsion,  but  they  are  also  firmly  fixed  in  position  to  reduce  damage  to  the  reduction  gear 
due  to  transmitted  vibrations.  Typically,  these  vibrations  are  considered  the  result  of  mass 
unbalance  in  the  machine's  rotating  parts,  however,  the  fluid  dynamics  of  the  air  flow 
within  the  engine  could  produce  forces  of  the  same  order  of  magnitude  as  those  forces 
caused  by  mass  unbalance.  As  a  result,  the  forced  vibrations  of  these  engines  must  be  as 
low  as  possible  to  reduce  premature  engine  failure. 

As  the  Navy  expanded  the  use  of  these  engines,  problems  developed  concerning 
residual  vibrations.  All  engines  delivered  to  the  Navy  met  specifications  for  vibration 
levels,  however,  some  engines  were  exceeding  vibration  levels  when  in  operation.  Gas 
generator  rotor  unbalance  was  identified  as  the  major  cause  of  higher  then  normal  vibration 
and  premature  wear.  As  a  result,  the  Trim  Balancing  Technique  was  developed  based  on 
the  Least  Squares  Influence  Coefficient  method  to  reduce  residual  unbalance.  The  great 
appeal  of  this  technique  is  that  the  whole  procedure  may  be  done  in  place  at  a  fraction  of  the 
cost  of  removal  and  replacement  of  the  engine.  This  technique  has  worked  very  well  and 
saved  the  Navy  millions  of  dollars  in  removal  and  replacement  costs,  but  it  is  not  100% 
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effective  nor  is  it  permanent.  With  increasing  operational  time  some  engines,  previously 
balanced,  show  an  increase  in  residual  vibration,  while  others  prove  incapable  of  being 
balanced  using  this  technique  [Ref.  1]. 

The  purpose  of  this  Thesis  is  to  examine  forces  contributing  to  these  vibrations 
other  then  simple  mass  unbalance.  Age  and  increased  operating  time  will  induce  ware 
and/or  shifting  of  parts  resulting  in  mass  fluctuations,  however,  it  is  possible  fluid  effects 
of  non-uniform  flow  through  the  engines  rotating  parts  could  produce  forces  approaching 
those  forces  caused  by  mass  unbalance.  Some  of  these  rotating  parts  that  could  cause  non 
uniform  flow  are  rotating  blades,  blade  tip  shrouds,  and  labyrinth  seals.  Such  forces  may 
even  be  correctable  by  current  mass  balancing  techniques  if  the  force  and  phase  of  the  fluid 
effects  are  compatible  with  the  modal  shapes  selected  for  balancing. 

This  thesis  will  develop  a  model  that  predicts  the  force  resulting  from  the  fluid 
effects  induced  by  a  non-circular  rotating  labyrinth  seal.  Many  of  these  seals  exist  in 
current  marine  engines,  and  possess  very  tight  sealing  gaps.  The  mean  radius  of  the  seal  is 
usually  very  large  compared  to  the  allowed  tolerance.  As  a  result,  even  a  small  defect  in 
seal  radius  compared  to  the  small  size  of  the  gap  will  generate  a  relatively  large  effect. 
Given  specific  operating  flow  conditions,  forces  can  be  predicted  that  will  rival  those 
allowed  by  acceptable  residual  mass  unbalance.  Given  various  geometry's  and  operating 
conditions  it  is  also  possible  to  predict  seal  flow  regimes. 


II.   BACKGROUND 

One  of  the  major  concerns  in  rotating  machinery  design  is  the  residual  unbalance 
due  to  the  mass  center  of  the  rotating  element  not  being  co-planer  with  its  axis  of  rotation. 
Such  a  situation  will  result  in  synchronous  lateral  vibration  on  the  rotor  bearing  elements. 
Under  certain  situations  the  magnitude  of  the  lateral  displacements  caused  by  this  vibration 
could  result  in  premature  failure  of  the  rotor  due  to  fatigue  or  catastrophic  failure  of  the  load 
bearing  elements.  Currently  mass  unbalance  is  thought  to  be  the  major  cause  of  such 
vibrations,  however,  bearing  wear,  non  uniform  thermal  expansion,  misalignment  between 
parts,  and  fluid  flow  within  the  machine  are  other  possible  causes  of  lateral  vibration. 
Much  work,  both  analytical  and  experimental,  has  been  done  to  detect  and  reduce  these 
unbalances.  In  the  1930's  dynamically  balancing  of  rotors  was  widely  practiced  on  test 
stands  with  moderate  results.  As  the  speed  of  the  machines  increased  this  method  proved 
less  effective.  Loewy  and  Piarulli  [Ref.  2]  point  this  out  and  discuss  the  advantages  of 
modal  balancing  and  self-balancing  as  effective  means  of  reducing  vibrations.  Modal 
balancing  works  well  with  the  lower  modes  as  long  as  the  harmonic  content  distribution  of 
the  imbalance  force  is  not  scattered  in  the  higher  modes.  Self-balancing  makes  use  of  two 
freely  swinging  weights  within  an  annular  track.  As  the  machine  rotates  the  weights  are 
free  to  take  equilibrium  azimuthal  positions.  At  speeds  below  critical  this  method  tends  to 
add  to  the  unbalance,  however,  above  critical  speed  the  weights  reduce  the  lateral  unbalance 
to  zero.  Both  these  methods  work  for  constant  speed  machines,  yet  they  are  not  practical 
for  variable  speed  machines  such  as  those  used  by  the  Navy. 

For  large  rotating  machinery  it  is  imperative  to  understand  the  natural  characteristics 
of  each  unit,  not  only  to  predict  critical  speeds  and  safe  operating  ranges,  but  also  be  able  to 
diagnose  danger  signs  in  machinery  history  to  avoid  catastrophic  failure.    Muster  and 


Stadelbaur  [Ref.  3]  discuss  the  necessity  of  dynamically  balancing  a  rotor  in  two  or  more 
axial  correction  planes  at  one  or  more  operating  speeds.  Den  Hertog  [Ref.  4]  discuses  the 
technique  of  shop  balancing  and  then  field  balancing  at  rated  speed.  Some  machines  at 
higher  critical  speeds  require  multi-planer  balancing  which  does  not  always  result  in  a 
balanced  rotor.  Loewy  and  Piarulli  [Ref.  2]  discuss  the  modal  balancing  theory  that  makes 
use  of  the  natural  frequencies  and  mode  shapes  of  the  shaft,  which  are  normally  needed  for 
super  critical  operations.  This  method  works  well  as  long  as  higher  harmonic  mode 
deflection  shapes  do  not  resemble  the  unbalance  deflection  shape. 

With  the  arrival  of  the  GE  LM-2500  Gas  Turbine  Engine,  the  United  States  Navy  has 
become  acutely  aware  of  the  difficulties  of  rotor  unbalance  and  the  resulting  decrease  in 
engine  service  life.  In  the  late  1970's  the  SPRUANCE  class  destroyers  and  the  OLIVER 
HAZARD  PERRY  class  frigates  began  fleet  operations  with  a  new  propulsion  plant 
centered  on  the  LM-2500.  Originally  this  engine  had  a  service  life  of  5000  hours,  yet  an 
alarming  rate  of  engines  could  not  meet  this  requirement  due  to  high  in  service  residual 
vibrations.  All  engines  met  required  vibration  levels  before  delivery,  therefore,  something 
was  causing  an  increase  in  vibrations  while  these  engines  were  in  operation.  As  a  result 
Thompson,  Badgley,  and  Raczkowski  [Ref.  1]  developed  a  trim  balancing  technique  that 
allows  for  the  balancing  of  the  engine  in  place  at  great  cost  savings.  Thompson  [Ref.  5] 
points  out  that  this  technique  can  not  always  be  used  effectively.  Some  engines  simply  are 
not  capable  of  being  balanced  whether  due  to  in  service  damage,  bowed  rotor,  or  some 
other  factor. 

One  such  factor  could  be  fluid  forces  acting  on  the  rotor  due  to  asymmetries  in  the 
rotor  airfoils  and  labyrinth  seals  rather  then  simple  mass  unbalance.  The  focus  of  this 
study  will  be  the  potential  contribution  to  synchronous  forces  from  non-circularity  in  the 
labyrinth  seals  and  what  the  magnitude  of  these  forces  could  be  given  various  design 
criteria.  With  the  tolerance  for  a  seal  radius  being  a  few  thousandths  of  an  inch  on  a  disk 


with  a  radius  of  one  to  two  feet  and  a  nominal  gap  clearance  on  the  order  of  thousandths  of 
an  inch,  there  exists  a  possibility  for  substantial  gap  variation.  It  is  the  mean  gap  clearance 
along  with  the  radius  tolerance  that  produces  the  gap  variation. 

A  labyrinth  seal  is  not  a  seal  at  all.  Given  a  proper  rotor/stator  configuration,  there 
is  no  contact  between  the  rotating  labyrinth  knives  and  the  stationary  casing.  The  propose 
of  the  labyrinth  is  to  create  a  path  for  the  fluid  to  travel  that  causes  a  substantial  change  in 
kinetic  energy  within  the  seal  gland  area.  This  change  in  kinetic  energy,  which  dissipates 
in  successive  chambers,  creates  areas  of  pressure  variation  upstream  of  the  gland  compared 
to  the  exit  pressure  of  the  gland,  which  seals  the  upstream  from  the  downstream 
conditions. 

A  perfect  operating  seal  with  no  intolerance  in  the  seal  gaps  has  uniform  pressure 
and  flow  momentum  changes  within  the  seal  gland  area.  There  is  no  resultant  radial  force 
on  the  shaft.  If  some  non  circular  intolerance  did  exist  on  the  seal  knife,  then  the  flow 
through  the  seal  would  not  be  uniform,  causing  areas  of  higher  mass  flow.  The  ultimate 
result  is  a  non  uniform  pressure  distribution.  If  a  non  uniform  pressure  distribution  exists 
in  the  gland,  it  is  a  simple  calculation  to  convert  it  to  a  radial  force  in  magnitude  and  phase 
relative  to  the  rotating  reference  system. 

Gas  turbine  rotors  are  typically  balanced  to  G2.5  limits  as  specified  by  ISO  1940 
[Ref.  6].  Given  these  constraints  the  goal  is  to  maintain  the  residual  mass  unbalance  force 
to  less  than  1  N  for  every  Kg  of  rotor  mass  at  10,000  RPM.  This  mass  unbalance  force  is 
described  as  F  =  mereco2,  where  m/e  is  the  equivalent  unbalance  (kg-m)  and  co  is  the 

shaft  rotational  speed  (rad/sec).  Once  the  fluid  force  is  known  it  is  compared  to  that  of  the 
mass  unbalance  force.  A  comparison  gives  some  indication  of  the  significance  such  fluid 
forces  could  have  on  rotor  unbalance  and  the  possible  necessity  to  consider  them  in  the 
design  process.  Under  certain  conditions  these  forces  may  present  no  significant  obstacle 
to  the  modal  balancing  techniques  currently  used.  If  the  force  is  of  significant  magnitude  or 


phase  to  impact  on  specific  mode  shapes,  correction  by  simple  adjustment  of  the  mass 
distribution  may  not  be  possible,  which  is  precisely  the  condition  experienced  in  some  LM- 
2500's  (Thompson  [Ref.  5]). 


III.     MODEL  FORMULATION 

A.       FLOW  EQUATIONS 

The  flow  in  a  labyrinth  seal  is  three  dimensional,  turbulent,  highly  unsteady,  and 
compressible.  Since  the  effects  of  this  complicated  flow,  not  the  details,  are  of  direct 
interest  for  this  problem,  a  simplified  lumped  parameter  fluid  model,  used  by  Kostyuk 
[Ref.  7],  Iwatsubo  [Ref.  8],  and  Millsaps  [Ref.  9],  will  be  introduced.  This  model  treats 
the  flow  in  the  labyrinth  chamber  as  a  one-dimensional  flow  in  the  azimuthal  direction,  but 
allows  for  axial  influx  and  efflux  due  to  the  pressure  driven  flow  through  the  seal  gaps. 
This  type  of  model  contains  the  dominant  physics  that  generate  a  non-uniform  pressure 
distribution  in  the  azimuthal  direction  within  the  gland.  The  mechanism  driving  this 
pressure  distribution  is  non-uniform  mass  and  azimuthal  momentum  exchange  between  the 
upstream  receiver  and  the  seal  gland.  Along  with  the  basic  model  simplifications,  several 
other  simplifications  are  made 

1 .  Ideal  gas.  Pv=RT 

2 .  Adiabatic  Process 

3 .  The  axial  flow  is  calculated  using  a  Bernoulli  equation  based  on  the  local  pressure 
difference  at  a  given  time  and  angular  location  using  the  average  density. 

4.  The  frictional  shear  stress  on  the  fluid  in  the  azimuthal  direction  within  the  gland 
can  be  calculated  using  the  Darcy  friction  law  for  1-D  flow. 

5 .  No  duct  acoustic  resonance  occur  in  the  range  of  interest. 

It  is  also  assumed  that  the  seal  is  fed  by  a  uniform  upstream  plenum.  This  condition  is  not 
strictly  met,  however,  the  perturbation  flow  will  induce  a  flow  redistribution  in  the 
upstream  volume,  helping  to  simplify  the  analysis.  Care  must  be  exercised,  however, 
since  under  certain  conditions  the  non  uniformity  in  the  upstream  flow  can  have  a  strong 


impact  on  the  flow  in  the  seal.  In  particular,  a  large  force  augmentation  is  possible  as  a 
result  of  the  upstream  non-uniformity.  An  analysis  of  this  phenomena  for  a  whirling  seal 
with  a  finite  upstream  volume  is  given  by  Millsaps  and  Martinez  [Ref.  10]. 

The  one-dimensional  continuity  equation  for  the  flow  along  the  azimuth  inside  the 
seal  gland  may  be  written  as 


Jt[pA]+i-k[pAV]+"2-"r0  e> 

s 


where  the  first  two  terms  are  the  normal  1-D  continuity  relation  for  a  deformable  duct.  The 
fluid  density  is  p ,  V  is  the  swirl  velocity  in  the  gland,  and  A  the  gland  cross  sectional  area 
normal  to  V.  The  remaining  terms  are  due  to  the  axial  mass  flux  over  the  outlet  and  inlet 
knives  respectively. 

The  one-dimensional  momentum  equation,  which  includes  the  rotor  and  stator 
frictional  shear  stresses,  Ts  and  xr  ,  is  similarly  written  as, 


d  r    ,„,       I     d 
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—  \pAV]  + pAVL    +  qnV-q,V.+  T  I- T  [l  +  2h)  + =  0  (2) 

dt  L  J     R    del  J     y2       Mi       s        ry  '     R    dd 

s  s 

The  inlet  swirl  is  given  by  Vt  ,  and  the  last  term  is  the  induced  pressure  gradient  force.  The 
total  incoming  mass  flow  is  m=pAs  w  where  As  is  the  effective  circumferential  gap  area 

and  w  is  the  jet  axial  velocity.  The  effective  inlet  area  seen  by  the  flow  is  the  real  flow  area 
multiplied  by  a  flow  coefficient  {fi ).  The  first  knife  coefficient  will  typically  be  fix  =  .65. 
The  resulting  mass  influx  per  unit  circumferential  length  ( q  , )  is  therefore, 

qx=llxpx8xVx  (3) 


A  similar  equation  applies  for  the  mass  efflux  per  unit  circumferential  length,  however,  the 
flow  coefficient  may  be  somewhat  higher  due  to  the  kinetic  energy  carry  over  effect. 

To  calculate  the  inflow  (q ,)  Bernoulli's  equation  is  used  with  the  density  being 
averaged  between  the  upstream  and  seal  gland.  This  yields, 


<?>  = 


iM, 


v*. 


\(p;-p2J  (4) 


where  P,  <5,   and  qx  are  functions  of  time  and  angle.  Similarly,  for  the  outlet  flow 


"^Jwr(p'-p^  m 


The  shear  forces  in  the  gland  from  both  rotor  and  stator,  which  are  responsible  for 
changing  the  average  swirl  of  the  fluid  in  the  gland,  are  calculated  by  using  the  Darcy 
friction  coefficients  law 


ts  =  ^Xssgn{V)lV2        and       rr  =  £  Xrsgn(coRx  -  V)(2h  +  l)(V-  (oRsf  (6) 
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The  friction  coefficients.  Av  and  Ar,  are  calculated  using  Blausius'  formula  with  the 
relevant  local  Reynolds  number,  based  on  the  gland  hydraulic  diameter  and  the  relative 
tangential  velocity  between  the  solid  surface  in  question  and  the  mean  swirl  velocity  within 
the  gland 


A     =0.3164<K^25)  (7) 


These  equations  will  be  used  to  determine  the  circumferential  flow  inside  a  spinning  seal 
that  does  not  have  perfect  circularity. 

B.         SEAL  IMPERFECTION 

A  labyrinth  seal  is  nominally  circular  unless  there  are  small  deviations,  which  could 
be  due  to  machining  tolerances  as  well  as  in  service  degradation.  A  very  small  "waviness" 
in  the  height  of  the  seal  knives  can  introduce  substantial  relative  variations  in  flow  area, 
especially  at  operating  conditions  where  seal  clearances  are  the  tightest. 

Letting  the  circumferentially  averaged  gap  for  the  first  knife  be  denoted  by  5* ,  and 
the  local  gap  be  8{(9),  where  9  is  the  angular  coordinate,  the  resulting  difference  (error)  is 
<?, .  The  local  gap  for  a  stationary  seal  can  be  expressed  as 

Sl(0)  =  8*l+el(9)  (8) 

The  net  radial  force  on  the  rotor  due  to  a  pressure  distribution  around  the  seal  gland 
is 
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Fr=-7TljP(9)e,ed9  (9) 


From  orthogonality,  only  the  first  harmonic  component  of  P{9)  will  contribute  to  the  radial 
force.  Since  the  perturbation  analysis  to  follow  is  purely  linear,  it  is  sufficient  to  retain  the 
first  harmonic  of  the  seal  geometric  imperfection  to  obtain  the  first  harmonic  in  pressure 
perturbation.  The  first  harmonic  of  the  axial  flow  area  is  given  by 
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e  =  -  \S(9)cos(9)d9  (10) 
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For  convenience,  the  extra  flow  area  can  be  concentrated  at  a  single  location  and 
would  appear  as  a  notch.  Mathematically  this  can  be  treated  as  a  Dirac  delta  function. 
Alternatively,  this  imperfection  is  equivalent  to  a  circular  seal  whose  center  is  displace  a 
distance  e  from  the  center  of  rotation  and  land  centroid.  Figure  1  shows  a  typical  non- 
uniform sealing  gap  for  a  single  knife,  along  with  the  equivalent  first  harmonic  and  notch 
representations.  Figure  2  shows  the  2-D  geometry  of  a  single  gland  seal  with  a  notch 
imperfection  modeled  as  a  single  notch.  Figure  3  shows  both  the  first  and  second  knives  in 
a  3-D  representation,  with  the  first  knife  imperfection  being  at  a  different  location  as  the 
second  knife  imperfection.  When  spinning,  the  first  knife  non-dimensional  gap  distribution 
is  given  by 


5,(0,/)  _ 


s; 


=  l  +  £,en°-lon  (11) 


Where  e,  is  the  amplitude  of  the  first  harmonic  deviation  from  circularity.  Similarly,  for  the 
second  knife  the  imperfection  is  represented  by 


^M  =  1  +  ^<*-«>  (12) 

82 


Finally,  a  representation  for  the  seal  gland  not  having  the  same  center  as  the  casing  is 
made.  The  height  of  the  gland  is  given  by 


4  =  1  + V""'  (13) 


1  1 


Note  that  e, ,  £2  and  e3  are  in  general  complex  so  that  they  contain  phase  information  on 
the  location  of  the  imperfections  as  well  as  the  amplitude  of  the  imperfection.  A  perfect 
casing  will  be  considered  here.  A  first  harmonic  deviation  in  the  land  would  produce  a 
steady  force  on  the  rotor  in  the  absolute  frame.  Since  the  magnitude  of  this  force  would  be 
small  compared  to  the  weight  of  the  rotor,  it  is  probably  of  little  importance  and  will  not  be 
considered  in  this  model. 

The  imperfection  generating  the  forces  does  not  have  to  be  due  to  sealing  gap  area 
variations.  Flow  coefficient  asymmetries  due  to  tip  imperfections  can  also  generate  mass 
flow  perturbations  and  hence  radial  forces.  For  the  purposes  of  this  model,  only  forces 
generated  from  seal  gap  non-circularity  will  be  considered. 

C.         SOLUTION  TECHNIQUE 

An  approximate  solution  for  the  azimuthal  pressure  in  the  seal  will  be  found  by  a 
linear  perturbation  analysis.  The  zeroth  order  solutions  for  the  case  of  a  perfect  seal,  which 
are  independent  of  time  and  angle,  are  found  from  Equations  (1)  and  (2)  by  noting 
d[]l  dt  =  0  and  d[]/  dd  =  0  .  All  zeroth  order  solutions  will  be  denoted  with  []*.  The 
steady  pressure,  density  and  flow  rate  are  extracted  from  the  zeroth  order  solution  to 
Equation  (1) 


P  = 
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Equation  (2)  is  used  to  solve  for  the  zeroth  order  swirl  velocity  (V*).  For  fixed  friction 
factors  the  resulting  equation  is  of  a  quadratic  form  in  V' .  Since  the  friction  factors  are 
functions  of  V\  in  that  they  are  a  function  of  the  local  Reynolds  number,  the  following 
quadratic  equation  in  V*  must  be  solved  by  iteration.  The  inlet  swirl  velocity  (  V()  is  used 

as  an  initial  guess 


■[Xj-Xr(l  +  2h)]v*2  + 
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V   - 


q'Vi  +  P-Ai.(K\  +  2h)(D2R~ 


=  0 


(17) 


With  the  flow  conditions  for  the  problem  of  the  perfect  rotor  obtained,  a 
perturbation  solution  for  the  imperfect  rotor  will  be  sought.  Since  it  is  assumed  that  the 
pressure  and  velocity  inside  the  gland  can  be  represented  as  the  zeroth  order  values  plus 
small  perturbations,  which  are  functions  of  (t,  9),  they  may  be  represented  as 


P  =  P*[l  +  4{d,t)] 
V  =  V'[\+ri(0,t)] 


(18) 


The  first  harmonics  of  these  being  functions  of  time  and  angle  can  be  expressed  as 
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When  these  expressions  are  substituted  into  Equation  (1)  and  linearized,  the  continuity 
equation  becomes 
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Similarly,  the  linearized  momentum  equation,  Equation  (2),  becomes 
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Equations  (21)  and  (22)  are  a  set  of  linear  complex  algebraic  equations,  easily 
solved  to  yield  the  complex  amplitudes  of  the  pressure  and  velocity  perturbations.  The 
pressure  amplitude  is  given  in  terms  of  the  real  and  imaginary  parts,   ^   and    £3, 

respectively  as 


£  =  £«+'£ 


(23) 


The  amplitude  of  the  pressure  perturbations  is 


i|i=Vs+fi 


(24) 
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and  the  phase  relative  to  the  reference  imperfection  is 


V. 


=  tan 


V 

{%*) 


(25) 


To  calculate  the  force  on  a  rotor  due  to  the  non-uniform  pressure  distribution, 
which  is  now  represented  by  the  amplitude  of  the  pressure  perturbation,  it  is  now  a  simple 
task  of  converting  this  pressure  to  force.  This  magnitude  is  all  that  is  sought  in  this 
analysis,  however,  it  should  be  noted  that  phase  information  can  also  be  extracted,  which 
may  prove  important  in  modal  vibration  analysis  and  balancing. 
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VI.     MODEL  PREDICTION 

A.  SEAL  FLOW  CONDITIONS 

The  solutions  to  Equations  (21)  and  (22)  depend  on  the  zeroth  order  parameters. 
These  conditions  depend  on  nominal  seal  geometry,  inlet  flow  conditions  ( P.  and  Vt),  and 

speed  of  the  rotor.  Figure  4,  which  shows  solutions  to  Equation  (17)  based  on  a  wide 
variation  in  nominal  seal  geometry,  inlet  flow  conditions  (P,  and  Vt),  and  speed  of  the 

rotor,  provides  non  dimensional  interrelationships  among  these  parameters  for  a  fixed  seal 
geometry  typical  of  modern  gas  turbines  (namely  the  sealing  height  and  pitch  being  equal). 
Figure  4  is  a  map  that  indicates  what  flow  "regime"  occurs  for  a  given  set  of  conditions. 
For  example,  if  the  axial  flow  rate  through  the  seal  is  small  and  the  difference  between  the 
peripheral  speed  of  the  rotor  (Q)RS)  and  the  inlet  swirl  velocity  (V.)  is  large,  then  the 

fractional  change  in  swirl  will  be  large.  On  the  other  hand,  if  the  flow  rate  is  large,  due  to  a 
large  pressure  gradient  or  large  sealing  gaps,  and  the  wheel  speed  is  about  the  same  as  the 
inlet  swirl,  there  will  be  a  relatively  small  fractional  swirl  change.  If  E  is  about  1.5,  there  is 
little  tendency  for  the  swirl  to  change,  as  would  be  expected,  since  this  is  the  "equilibrium" 
swirl  value  that  would  be  obtained  in  the  absence  of  an  axial  through  flow. 

B.  RADIAL  FORCE  PREDICTIONS 

Results  for  particular  cases  will  now  be  presented  to  illustrate  the  use  of  this  model 
in  the  prediction  of  synchronous  radial  rotor  forces.  Table  1  shows  the  seal  geometry  used. 
The  seal's  operating  conditions  and  geometry  are  typical  for  an  LM-2500  gas  generator. 
The  dominant  parameters  controlling  the  forces  are: 
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1 .  The  pressure  difference  across  the  seal  knives. 

2.  The  inlet  swirl  velocity. 

3 .  The  swirl  change  through  the  seal. 

4.  The  magnitude  of  the  seal  imperfection  relative  to  the  nominal  gap. 

5 .  The  angle  of  the  second  seal  imperfection  relative  to  the  first. 

The  general  scaling  behavior  for  the  radial  synchronous  force  is  approximately  proportional 
to  the  inlet  total  to  exit  static  pressure  difference  (Pi  -P„)  as  would  be  expected  for 

incompressible  flow.  Some  nominal  deviation  from  this  scaling  occurs  at  higher  pressure 
ratios  due  to  compressibility  effects.  Due  to  the  linearity  of  the  model,  the  predicted  forces 
will  scale  directly  with  the  magnitude  of  the  imperfection  used  for  the  seal.  When  the 
magnitude  of  the  seal  non-circularity  doubles  for  a  given  nominal  seal  gap,  the  force 
doubles  accordingly.  Similarly,  for  a  fixed  non-uniformity,  the  force  scales  inversely  with 
the  nominal  seal  gap.  With  the  force  scalings  proportional  to  the  geometry,  it  will  be 
possible  to  non-dimensionalize  and  present  the  force  as 


F=^— =  ^^  (26) 

APlRs£t      APlRse} 


1.  Single  Knife  Imperfection 

The  magnitude  of  the  radial  force  due  to  an  imperfection  in  the  first  sealing 
knife  is  shown  in  Figure  5.  The  force  coefficient  is  shown  verses  rotational  frequency  for 
several  values  of  inlet  swirl  velocities.  As  the  inlet  swirl  increases,  the  frequency  at  which 
the  minimum  force  occurs  also  increases.  To  obtain  a  better  understanding  of  the 
phenomenon  taking  place,  Figure  6  is  shown  with  the  exact  conditions  as  those  for  Figure 
5,  except  the  effects  of  friction  are  taken  out  of  the  model.  With  no  friction  there  is  no 
mechanism  for  changing  the  swirl  of  the  incoming  flow  and  as  a  result  the  minimum  force 
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occurs  when  the  rotational  speed  of  the  rotor  equals  the  swirl  speed  in  the  gland.  This 
suggests  that  the  minimum  force  occurs  when  the  rotor  speed  is  equal  to  the  swirl  speed  in 
the  gland.  Returning  to  Figure  5,  the  effects  of  friction  are  not  only  impacting  on  the  swirl 
speed  in  the  gland,  but  the  momentum  change  of  the  flow  absorbs  some  of  the  energy  from 
the  force  that  could  have  been  imparted  to  the  rotor.  At  the  minimum  force  condition  the 
gland  swirl  speed  is  the  same  as  the  rotor  speed,  meaning  there  is  no  relative  motion 
between  them.  The  friction  within  the  gland  always  changes  the  speed  of  the  incoming 
swirl  depending  on  the  relative  velocity  between  the  inlet  swirl  and  the  rotor  speed.  Figure 
7  shows  the  phase  or  relative  position  where  the  force  acts  on  the  rotor.  The  phase  shift 
corresponds  to  that  frequency  where  the  force  is  a  minimum.  The  phase  is  a  function  of 
rotor  speed. 

2.  Duel  Knife  Imperfection 

In  this  case  the  same  conditions  exist  as  in  Table  1,  however,  both  the  front 
knife  and  the  rear  knife  are  capable  of  having  a  non-circularity.  These  non  circularities  are 
also  capable  of  having  a  phase  angle  between  them,  which  may  very  from  0  to  180 
degrees.  Figure  8  shows  the  effect  of  the  phase  angle  of  the  back  knife  notch  relative  to  the 
front  notch  as  a  function  of  rotor  speed.  The  size  of  each  imperfection  is  the  same,  and  the 
inlet  swirl  is  maintained  at  150  m/s.  Results  are  shown  for  the  angle  between  the 
imperfections  varying  from  0  to  180  degrees.  The  force  tends  to  be  less  when  the  notches 
are  at  the  same  location.  This  is  due  to  a  local  increase  in  pressure  at  the  front  seal 
imperfection  being  partially  negated  by  a  local  decrease  in  pressure  at  the  rear  seal 
imperfection.  It  should  also  be  noted  that  the  maximum  force  occurs  when  the  rotor  speed 
equals  the  speed  of  the  gland  swirl,  which  is  not  consistent  with  the  trend  of  having  a 
minimum  force  at  this  point.  This  case  can  be  described  as  a  bowed  rotor.  Even  though  the 
pressure  effects  tend  to  cancel  they  are  still  significant.   When  the  imperfections  are  not 

18 


lined  up  the  forces  are  larger  implying  the  pressure  effects  of  the  front  and  back 
imperfections  are  additive. 

C.         MULTIPLE  SEALS 

Up  to  this  point,  the  forces  from  single  gland  labyrinth  seals  have  been  considered. 
The  analysis  that  has  been  presented  can  be  readily  extended  to  the  case  where  there  are 
multiple  sealing  knives  within  a  single  labyrinth  seal.  The  results  for  multi-cavity  seals 
show  similar  trends  to  those  of  a  single  gland  seal.  For  a  single  gland  labyrinth  seal  it  has 
been  shown  that  two  complex  algebraic  equations  will  predict  the  fluid  forces  within  the 
gland.  This  analysis  can  be  easily  expanded  for  N  number  of  chambers  within  a  multi- 
chamber  labyrinth  seal.  The  extension  of  this  analysis  will  yield  2-N  equations  with  2-N 
unknowns.  This  is  not  treated  explicitly  here. 

Typically,  a  gas  turbine  possesses  dozens  of  locations  in  the  secondary  flow  path 
where  labyrinth  seals  are  placed.  In  this  case  there  are  multiple  seal  glands  which  are  not 
coupled  along  the  flow  path.  It  is  important  for  a  designer  to  be  able  to  incorporate  the 
radial  forces  arising  from  imperfect  seals  at  different  axial  locations.  The  question  that 
arises  is:  How  does  one  "add"  the  forces?  Two  special  cases  will  now  be  discussed. 

1.  In  Phase 

If  the  imperfections  of  all  the  rotor  seals  have  the  same  angular  location  then 
the  forces  are  added.  This  would  be  the  case  if  the  rotor  were  slightly  bowed,  or  if  the  seals 
were  all  machined  together  with  a  slight  off-set  on  the  lathe.  The  result  being  that  all  the 
forces  would  occur  at  the  same  location.  This  can  be  expressed  as 

N 

^=XF-  =  yv/r  (27) 

(=1 


19 


Where  FT  is  the  total  radial  force,  Ft  is  the  force  from  the  ith  seal  and  N  is  the  number  of 
seals.  Using  the  conditions  of  Table  1  and  assuming  16  stages  with  a  uniform  non- 
circularity  of  1,  the  model  predicts  the  resulting  force  will  be  515  N.  A  non-circularity  of 
one  is  very  high  for  a  modern  gas  turbine,  but  given  the  tight  tolerances  it  is  not 
impossible.  A  non-circularity  of  .1  is  quite  possible  under  normal  operating  conditions  and 
would  yield  a  force  of  5 1 .5  N  on  the  rotor. 

2.  Random  Phase 

If  the  imperfections  are  random,  then  there  is  no  phase  correlation  between 
the  imperfection  at  one  location  and  any  other  location.  It  may  appear  in  this  case  that  the 
forces  may  totally  cancel  but,  this  is  analogous  to  the  case  of  a  random  walk.  It  can  be 
shown  that  the  most  probable  magnitude  of  such  an  addition  of  forces  with  random  phases 
would  scale  with  the  square  root  of  the  number  of  seals. 


FT  =  F^N  (28) 

Given  the  same  conditions  as  the  "bowed  rotor"  case,  the  model  predicts  the  force  on  the 
rotor  for  a  non-circularity  of  1.0  will  be  128  N.  Similarly,  if  the  non-circularity  were  .1  the 
force  on  the  rotor  would  be  12.8  N.  It  is  expected  that  the  net  radial  force  will  be 
somewhere  between  these  limits  of  completely  in  phase  and  totally  random. 
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D.         FLUID  FORCE  COMPARED  TO  UNBALANCE  FORCE 

The  fluid  force  must  be  of  the  same  order  of  magnitude  as  the  force  due  to  the 
residual  mass  unbalance  if  it  is  to  be  considered  important  in  the  dynamics  of  the  rotor. 
From  ISO  1940  [Ref.  6]  the  acceptable  radial  force  which  corresponds  to  a  G2.5  balancing 
grade  is  roughly  0.1  Newtons  per  kilogram  of  rotor  mass  at  1000  RPM  and  about  1.0 
Newtons  per  kilogram  at  10,000  RPM. 

Since  the  LM-2500  engine  has  been  used  to  model  the  geometry  used,  it  is 
appropriate  to  use  the  mass  of  the  rotor  to  make  comparisons.  The  mass  of  the  rotor  is 
about  500  kilograms.  The  normal  operating  speed  or  the  rotor  is  about  10,000  RPM.  With 
these  conditions  the  limiting  mass  unbalance  is  1.0  Newtons  per  kilogram.  In  the  cases 
previously  discussed  the  forces  predicted  convert  to  mass  unbalance  forces  of  1 .03  -  .3 
Newtons  per  kilogram  for  a  non-circularity  of  1,  and  0.1  -  0.03  Newtons  per  kilogram  for 
a  non-circularity  of  0.1.  These  predicted  unbalance  forces  are  on  the  same  order  of 
magnitude  as  the  limiting  mass  unbalance  force. 

TABLE  1.  STANDARD  CASE  CONDITIONS 

DIMENSIONS  ATMOSPHERICS 

Rs  =  .5m  Pi  =110000  Pa 

L  =  .01m  Po=  100000  Pa 

h  =  .01m  T  =  295K 
<S,  =  .0001  m 
S2  =  .0001  m 


CO  =  1000  rad/sec 
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V.     CONCLUSIONS 

A  lumped  parameter  fluid  model  was  presented  that  predicts  the  generation  of 
synchronous  radial  forces  due  to  a  rotating  imperfect  labyrinth  seal.  The  predictions  may 
be  based  on  knife  imperfections  for  the  front  knife,  the  back  knife,  or  both  with  any  angle 
between  them.  The  model  is  capable  of  predicting  the  magnitude  of  the  resulting  force  and 
the  angular  location  on  the  rotor.  The  flow  characteristics  produced  by  the  model  are 
consistent  with  those  expected,  in  that,  the  non-dimensional  flow  regimes  are  consistent 
with  physical  behavior. 

The  forces  predicted  by  this  model  are  of  the  same  order  of  magnitude  as  those 
recommended  for  some  common  balancing  grades.  The  phase  predicted  could  also  have  an 
impact  on  the  modal  balancing  techniques  currently  used,  since  the  phase  of  these  forces  is 
a  function  of  rotor  speed  not  rotor  mass.  This  behavior  must  be  incorporated  into  the  trim 
balancing  technique  if  the  balancing  masses  are  to  account  for  the  residual  fluid  forces. 

The  model  predicts  that  a  minimum  force  condition  is  always  consistent  with  a  rotor 
speed  that  matches  the  gland  swirl  speed.  This  implies  a  design  criterion  calling  for  an  inlet 
swirl  speed  that  matches  rotor  speed.  Such  a  condition  will  reduce  residual  fluid  forces 
imparted  to  the  rotor.  Currently  swirl  breaks  are  used  in  engine  design  to  reduce  inlet  swirl 
and  help  reduce  rotor  instability.  It  is  possible  that  while  swirl  breaks  help  solve  the  rotor 
instability  problem,  they  may  add  to  the  synchronous,  forced  vibration  problem.  This 
condition  is  not  yet  manifesting  itself  since  rotor  mass  is  still  relatively  high.  With 
increasing  advances  in  materials  and  rotor  bearing  design,  rotor  mass  will  decrease  while 
machine  power  density  increases.  At  some  point  the  fluid  forces  produced  by  component 
non  circularity  could  become  more  significant  than  the  mass  unbalance  forces.  Under  such 
conditions  it  may  not  be  possible  to  balance  these  forces  with  mass  corrections,  meaning 
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proper  design  and  the  ability  to  predict  such  fluid  forces  based  on  given  tolerances  will  be 
imperative  to  rotor  balance. 
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VI.     RECOMMENDATIONS 

Many  simplifying  assumptions  have  been  made  in  order  to  reduce  this  highly 
complex  flow  pattern  down  to  a  model  that  could  make  some  predictions  consistent  with 
the  dominant  physics  taking  place  within  the  gland.  The  results  are  promising,  however, 
they  must  be  verified  with  experimental  data.  It  is  recommended  that  a  test  facility  be 
constructed  so  that  the  conditions  described  here  may  be  reproduced  on  a  test  platform.  It  is 
anticipated  that  discrepancies  will  arise,  however,  if  scaling  factors  can  be  applied  to  make 
the  model  accurate  then  it  will  be  possible  to  predict  what  design  tolerances  will  produce  in 
the  way  of  fluid  forces  on  the  rotor. 

The  phase  component  of  the  force  produced  was  briefly  touched  on  in  this  study, 
yet  it  may  prove  to  have  a  profound  effect  on  current  balancing  techniques,  since  the  phase 
of  this  force  is  a  function  of  rotor  speed  not  rotor  mass.  It  is  recommended  that  this  be 
studied  further. 

On  the  modern  gas  turbine  engine,  labyrinth  seals  are  not  the  only  source  of  fluid 
induces  synchronous  forces.  Each  rotating  blade  tip  is  structured  to  act  as  a  seal  in  order  to 
reduce  leakage,  hence  increasing  component  efficiencies.  Any  non  circularity  in  blade  tip 
radius  will  cause  non  uniform  radial  clearances  resulting  in  non  uniform  radial  forces.  As 
this  model  has  shown,  blade  to  blade  variation  in  radial  forces  could  affect  rotor  balance. 
These  radial  forces  at  the  blade  tips  should  be  studied  for  possible  contributions  to 
synchronous  force  vibrations  and  the  balancing  problem. 
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APPENDIX 


a) 


b) 


O 


Figure  1.    Three  equivalent  sealing- gap  representations  for  a  non-uniform 
single  labyrinth  knife:  a)  Exaggerated  non-uniformity, 
b)  Equivalent  first  harmonic,  c)  Equivalent  notch 


Figure  2.  Two  dimensional  view  of  notch  reference  frame  for  single  knife 
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Figure  3.  Three  dimensional  view  of  two  knife  labyrinth  seal 
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Figure  4.    Plot  for  nondimensional  interrelationship  for  seal  through  flow 
(Q,  shaft  spinning  speed  (E),  and  swirl  change  (D)  in  the 
labyrinth  gland. 
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